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SYSTEM AND METHOD FOR IDENTIFYING OPTIMIZED BLOOD SIGNAL IN 
MEDICAL IMAGES TO ELIMINATE FLOW ARTIFACTS 
Field of the Invention 

The present invention is directed to the identification of blood signals in medical 
5 images, such as magnetic resonance (MR) images, and more particularly to such 
identification which overcomes the corruption of such signals by arterial flow artifacts. 
Description of Related Art 

Dynamic contrast enhanced MRI (dceMRI) has demonstrated considerable utility in 
both diagnosing and evaluating the progression and response to treatment of malignant 
10 tumors. dceMRI involves the periodic acquisition of Tl -weighted images before, during and 
after injection of a gadolinium labeled tracer such as gadopentetate dimeglumine. The change 
over time in signal intensity in a voxel or region of interest in this time series can then be 
related to tracer concentration. By making use of a two-compartment model, with one 
compartment representing blood plasma and the other abnormal extra-vascular extra-cellular 
15 space (EES), the observed uptake curves in tissue and plasma can be used to estimate various 
physiological parameters. 

The parameter of primary interest in this work is the volume transfer constant 
between blood plasma and EES, commonly referred to as A*"™ 5 . This parameter is related to 
both blood flow and endothelial permeability, and is therefore a good endpoint for estimating 
20 the blood supply available to a target malignancy. 

One of the primary challenges in estimating perfusion parameters is identifying an 
accurate plasma uptake curve. Using a theoretical curve ignores differences in injection rate 
and cardiac output, which can greatly reduce reproducibility. However, the MR signal in 
arteries is frequently corrupted by flow artifacts, with the result that regions of interest at 
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different points in the same artery or in other nearby vessels can provide grossly different 
uptake curves. 

This problem is illustrated in Figs. 1 A and IB. In Fig. 1 A, one section from subject 1, 
time 1, has had two small arteries identified therein as regions of interest (ROFs). The ROI's 
5 identified on the right and the left are identified as vessel 1 and vessel 2, respectively. Fig. 
IB shows the raw uptake curves for vessels 1 and 2, along with a calculated ideal arterial 
uptake curve calculated in accordance with the present invention (to be described below). 

Note that in Fig. IB, vessel 1 and vessel 2 show significantly different uptake curves. , 
Moreover, neither vessel shows the distinctive sharp peak and subsequent plateau 
10 characteristic of plasma enhancement following injection of a tracer bolus. 
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Summary of the Invention 

It is an object of the invention to allow the calculation of blood uptake curves with 
improved accuracy. It is another object of the invention to allow such calculation such that 
flow artifacts are eliminated. 
5 To achieve the above and other objects, the present invention is directed to a system 

and method for the identification of an optimized plasma signal which is intended to 
eliminate this source of measurement variability and thereby increase the sensitivity to 
change of perfusion parameter measurements. The uptake curve according to the present 
invention shows a greater enhancement peak than that of either vessel, a smoother plateau, 

10 and a more characteristic shape. 

The present invention allows the calculation of perfusion parameters in imaging 
technologies such as dynamic contrast enhanced MRI. Uptake curves are calculated for both 
tumor tissue and plasma. Inter-operator variability in the derived rate constant between 
plasma and extra-cellular extra-vascular space is assessed using semi-automated tumor 

15 margin identification with both manual and automated plasma identification. In addition, an 
assessment is made of the contribution to total variability made by differences in tumor 
margin identification and differences in plasma identification. Experimental results show a 
mean coefficient of variability (CV) for parameter measurement with manual plasma 
identification of 20.1%, with a mean CV for parameter measurement with automated plasma 

20 identification of 6.7%. Analysis shows that 67% of the variability in parameter measurement 
with manual plasma identification is attributable to differences in identified plasma signal, 
with the remainder attributable to differences in identified tumor margins. 

Manual plasma identification for perfusion parameter calculation is currently standard 
practice for both clinical and experimental purposes. Increased accuracy and sensitivity to 
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change can be achieved by making use of an automated method for plasma identification, as 
achieved by the present invention. 

In some cases it may be possible to identify automatically the tumor and the artery, 
for example by use of special pulse sequences, pre-contrast, which create a known contrast 
5 with respect to surrounding tissue. 

The present invention can be applied to a variety of imaging technologies, such as 
dynamic contrast enhanced CT. 

The present invention bears some similarities to the algorithm presented by Rijpkama, 
M, Johannes, H., et al, "Method for quantitative mapping of dynamic MRI contrast agent 
10 uptake in human tumors," JMRI, pp. 457-463, 2001. However, the present invention 
provides automation. Also, because the method of Rijpkema et al selects most or all of the 
arterial voxels, it is vulnerable to the arterial flow artifacts which the present invention is 
intended to eliminate. 
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Brief Description of the Drawings 

A preferred embodiment of the present invention will be disclosed in detail with 
reference to the drawings, in which: 

Fig. 1 A is a raw image with two arteries identified therein; 
5 Fig. IB is a plot of blood uptake curves from the arteries identified in Fig. 1 A, as well 

as an optimal arterial uptake curve calculated in accordance with the preferred embodiment; 

Fig. 2 is a flow chart showing the process carried out in the preferred embodiment; 

Fig. 3 is a sample image from a perfusion data set used to test the preferred 
embodiment; 

10 Fig. 4 is a scatterplot of lC rans values using manual and automatic plasma 

identification; 

Fig. 5 is a trend line using manual plasma identification for a subject; 
Fig. 6 is a trend line using automated plasma identification for the same subject; and 
Fig. 7 is a block diagram of a system on which the preferred embodiment can be 
15 implemented. 
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Detailed Description of the Preferred Embodiment 

A preferred embodiment of the present invention, as well as experimental results, will 
now be set forth in detail with reference to the drawings. 

The process according to the preferred embodiment will be disclosed with reference 
5 to the flow chart of Fig. 2. Tumor margins are identified in step 202 using Geometrically 
Constrained Region Growth (GEORG). This technique requires a user to place a seed or 
string of seeds within each desired structure throughout the volume using one or more mouse 
clicks. The seed regions then expand into neighboring voxels provided that two constraints 
are satisfied: the grayscale value of the neighboring voxel must have a high probability of 
10 falling within the statistical distribution defined by all currently included voxels, and 
inclusion of the neighboring voxel must not cause the shape of the included region to deviate 
excessively from the a priori regional shape model. Once initiated, the expansion process 
continues until a stable boundary has been established. The resulting contour is then 
converted into a snake which can be interactively corrected by the analyst if the initial result 
15 is sub-optimal. 

Alternatively, step 202 could be performed by automatic identification of the tumor 
and the artery. One technique for doing so is by the use of special pulse sequences, pre- 
contrast, which create known contrast with respect to the surrounding tissue. 

After step 202, an optional step 203 of motion correction between time frames may be 
20 performed. Step 203 may be used to correct for breathing and other patient motions. 
Techniques for inter-frame motion correction are known in the art and will therefore not be 
disclosed in detail here, although their use in the context of the present invention is novel. 

After identifying the tumor margins, the analyst identifies a region of plasma in step 
204, preferably in an artery in close proximity to the tumor. This is done using manual tracing 
25 with a computer mouse. An alternative way to carry out step 204 is by using knowledge of 
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the timing of the injection and the start of scanning rather than the user-identified plasma 
region. Similarly to step 202, step 204 may be followed by an optional step 205 of inter- 
frame motion correction. 

At this point, the identified plasma region is used for parameter calculation, as 
5 described below. In addition, the identified plasma region is used in to initialize an automated 
search algorithm whose intent is to identify an optimized plasma signal for the data set under 
consideration. Each voxel in the data set is assigned a score in step 206 based on time point 
of maximum uptake, slope at maximum uptake, peak value, and conformance to a gamma 
variate curve. The highest scoring twenty-five voxels in the data set are then assigned in step 

10 208 to the ideal plasma region of interest. Thus, an optimized plasma signal is derived, the 
signal being optimized to eliminate flow artifacts. 

After plasma has been identified by either manual or automated means, uptake curves 
are generated for both tumor and plasma in steps 210 and 212, respectively. These are 
designated C t (t) and C p (t), respectively. In the interests of noise reduction, both plasma and 

15 tumor data are fit to gamma variate curves. The vascular bed is modeled as a linear system, 
such that: 

C,(0 = C,(0**(0 (0 
with impulse response h(i) given by: 

hify^K^e'^ (2) 
20 where k ep is the rate constant between the EES and blood plasma. Given C,(/) and C p (t), K! rans 
and k ep are estimated in step 214 using a gradient-descent energy minimization scheme. Local 
minima are avoided through the use of multiple instantiations with different initial parameter 
settings. In step 216, the I?™™ values can be mapped or imaged, or they can be superimposed 
on one of the original scan planes, so that the user can view the higher and lower K! rans values 
25 within the tumor. 
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Experimental results will now be presented. The experiments involved in this study 
were intended to assess the reproducibility of perfusion measurements using manual and 
automated plasma identification, and to determine the percentage of measurement variability 
due to differences in tumor margin and plasma region of interest, respectively. Experimental 
5 data were derived from three dogs with naturally occurring mammary tumors. Each animal 
was imaged three times over a period of 12 weeks. Images for this study were acquired using 
a GE 1 .5T LX/CV scanner. Three slices through each tumor were acquired using a cardiac 
coil. Perfusion images used a GRE pulse sequence with a repetition time of 20ms, echo time 
of 1ms, and a flip angle of 40 degrees. Imaging time for each image set was seven seconds, 
10 with a two second scanner delay, yielding temporal resolution for the data set of nine 
seconds. The reconstruction matrix was 256x192, FOV was 140mm, and slice thickness was 
4mm. 

Fig. 3 shows a sample image from the perfusion data set for Subject 1 at Time 3. The 
circles on either side of the tumor are phantoms. These were worn on a belt around the animal 

1 5 during imaging and were used for motion estimation and correction. 

Because a primary aim of the testing of the preferred embodiment was the assessment 
of inter-operator variability, four analysts were trained in the use of the analysis software. All 
analysts were also trained in the appearance of canine mammary tumors and the selection of 
appropriate plasma regions using images from animals not included in this study. Each 

20 analyst was then asked to identify and delineate both tumor and plasma in each of the nine 
included data sets. When identifying plasma, the analysts had the option to view the uptake 
curve for the currently selected region at any time, and to erase, modify or replace the 
currently selected region. In this way each analyst was able to manually select a reasonably 
optimized plasma region. 
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Once all regions of interest were delineated, K! rans values were calculated first using 
the regions of interest as identified by the analyst, and then using the analyst's tumor 
identification with the automatically identified plasma uptake curve. By comparing the 
variance seen between analysts using manually identified plasma with that seen between 
5 analysts using the automatically identified plasma, which was identical across analysts, it was 
possible to isolate variability related to plasma signal from that related to differences in tumor 
margin identification. 

Coefficients of variability in measurement of K* rans among the four analysts, defined 
as measurement standard deviation divided by measurement mean, were calculated separately 

10 for manual and automatic plasma identification, and for each of the nine cases examined. For 
the nine manual plasma identifications, coefficients of variability ranged from 3.1%. to 
39.2%, with a mean of 20.1% and a median value of 21.5%. For the nine automated plasma 
identifications, coefficients of variability ranged from 3.1% to 11.8%, with a mean of 6.7% 
and a median value of 6.2%. 

15 Bearing in mind that the same tumor margins were used for both the automated and 

the manual plasma calculations, it can be generally surmised that approximately two-thirds of 
the variability seen in the manual measurements was a result of differing plasma signal 
identifications, with the remaining one-third attributable to differing tumor margin 
identifications. It should be noted that the variability attributed to differences in tumor 

20 margins is similar to that reported previously for volume measurements of lung tumors using 
GEORG. 

An examination of a scatterplot of K* rans measurements using manual vs. automatic 
plasma identification, as shown in Fig. 4, shows that the correlation between the two 
measures is reasonable given the high variability of the manual measurements. It also shows 
25 a slope of 0.874, indicating that on average the manual measurement gives a somewhat 
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higher estimation of A^ ranj than the automatic measurement. This is as expected, since the 
general effect of flow artifacts will be to reduce the apparent plasma enhancement, thereby 
exaggerating the proportion of tracer apparently passing into the EES. 

Fig. 5 shows a trend line of I? rans using manual plasma identification for subject 1 . 
5 The subject cannot be confidently classified as declining until time 3. 

Fig. 6 shows a trend line oiK! rans using automated plasma identifcation for subject 1. 
The subject can be confidently classified as declining after time 2. 

An examination of the trend over time in f? rans for subject 1 using manual and 
automated plasma identification, as shown in Figs. 5 and 6, highlights the value of the 

10 reduced measurement variability afforded by the automated process. Although both trend 
lines indicate that vascular perfusion for this tumor is declining over time, higher variability 
makes that assumption statistically insupportable at time two for manual plasma 
identification, and marginally supportable after time three. Using automated plasma 
identification, however, this subject may be confidently classified as declining after time 2. 

15 The difficulty of identifying a suitable plasma signal is typically greater in smaller 

animals such as the dogs used in study than in humans. This is due to small animals' higher 
blood velocity, which exaggerates flow artifacts in the arteries, as well as to the lower signal 
to noise ratio that is achievable when imaging smaller anatomy. The values given above for 
parameter variability due to differences in plasma identification should be considered an 

20 upper limit when estimating likely variability in human studies. 

The preferred embodiment, as well as other embodiments, can be implemented on 
hardware such as the system 700 shown in Fig. 7. An input 702 receives image data, which 
can be previously stored, received from an imaging device, or communicated from a remote 
location. A CPU 704 receives the imaging data from the input and performs the operations 

25 disclosed above. The CPU interfaces with the analyst through a display 706, a keyboard 708 
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and a mouse or other pointing device 710; the utility of the mouse 710 for identifying seeds 
and regions has been disclosed above. The results can be output to any suitable output 712, 
e.g., a printer, a storage device, or a communication device for communicating the results to a 
remote location. 

5 While a preferred embodiment of the present invention has been set forth above, those 

skilled in the art who have reviewed the present disclosure will readily appreciate that other 
embodiments can be realized within the scope of the present invention. For example, the 
utility of the present invention is not limited to MRI. Also, any suitable hardware could be 
used; for example, instead of a mouse, a light pen, stylus and pressure pad, or other suitable 
10 pointing device could be used. Furthermore, any suitable technique for identifying tumor 
margins could be used. Therefore, the present invention should be construed as limited only 
by the appended claims. 
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